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ABSTRACT 

It has long been known how to analytically relate the clustering properties of the collapsed structures 
(halos) to those of the underlying dark matter distribution for Gaussian initial conditions. Here we 
apply the same approach to physically motivated non-Gaussian models. The techniques we use were 
developed in the 1980s to deal with the clustering of peaks of non-Gaussian density fields. The 
description of the clustering of halos for non-Gaussian initial conditions has recently received renewed 
interest, motivated by the forthcoming large galaxy and cluster surveys. For inflationary-motivated 
non-Gaussianites, we find an analytic expression for the halo bias as a function of scale, mass and 
redshift, employing only the approximations of high-peaks and large separations. 

Subject headings: cosmology: theory, large-scale structure of universe - galaxies: clusters: general - 
galaxies: halos 



1. INTRODUCTION 

Constraining primordial non-Gaussianity offers a pow- 
erful test of the generation mechanism of cosmologi- 
cal perturbations in the early universe. While stan- 
dard single-field models of slow-roll inflation lead to 
small departures from Gaussianity, non-standard sce- 
narios allow f or a larger level of non-Gaussianity 
(|Bartolo et al.l (120041 ) and references therein). The 
standard observables to constrain non-Gaussianity are 
the cosmic microwave background and large-scale struc- 
ture. A powerful technique is based on the abundance 
jMatarrese et al.' 2000; Verde ct al. 2001; LoVcrde et al.' 
[2007; Robinson & Baker 2000; Robinson et al. 200a) 
and clustering (iGrinstein fc Wiselll986l : iMatarrese et al.l 
119861 : iLucchin et al.lll988f) of rare events such as dark 
matter density peaks as they trace the tail of the 
underlying distribution. These theoretical predictions 
have been tested against numerica, l N-body simula- 
tions (iKang et all |2007t jGrossi et al.l l2007t iDalal et all 
120071 1. 103,^1 et all ()2007n showed that primordial non- 



Gaussianity affects the clustering of dark matter halos 
inducing a scale-dependent bias. This effect will be use- 
ful for constraining non-Gaussianity from future surveys 
which will provide a large sample of galaxy clusters over 
a volume comparable to the horizon size (e.g., DES, 
PanSTARRS, PAU, LSST, DUNE, ADEPT, SPACE, 
DUO) or mass-selected large clusters samples via the 
Sunyaev-Zel'dovich effect (e.g., ACT, SPT), considered 
alone or via cross-correlation techniques (e.g., ISW, lens- 
ing). 

Here, we resort to results and t echniques develop ed in 
the 1980s (IGrinstein fc Wisdll986l : [Matarrese eral1ll986l : 
ILucchin et al.lll988D to extend this work and derive an ac- 
curate analytical expression for halo bias, in the presence 
of general non-Gaussian initial conditions, accounting for 
its scale, mass and redshift dependence. 

^ Dipartimcnto di Fisica "G. Galilei", Universita degli Studi 
di Padova and INFN, Sezione di Padova, via Marzolo 8, 35131, 
Padova, Italy; sabino.matarrese@pd.infn.it 

2 ICREA and Institute of Space Sciences (GSIC-IEEC), UAB, 
Barcelona 08193, Spain 

^ Department of Astrophysical Sciences, Peyton Hall, Princeton 
University, Ivy Lane, NJ 08544, USA ; verde@ieec.uab.es 



2. HALOS AS PEAKS OF THE DENSITY FIELD 

Halo clustering is generally studied by assuming that 
halos correspond to regions where the (smoothed) linear 
density field exceeds a suitable threshold. This amounts 
to modeling the local halo number density as a theta 
(step) function 

Ph,R(x, 2/) = [(5fl(x, z/)- A,]=0 [5n{-^)-5c[zf)] , (1) 

modulo a constant normalization factor which is irrele- 
vant for the calculation of correlations. Here R denotes 
a smoothing radius which defines the halo mass M by 
M = n„,fi3H§/{8TTG){A/3)TrR^, with n,n,Q denoting the 
present-day matter density parameter, Hq the present- 
day Hubble parameter and G Newton's constant. The 
threshold Ac is the linearly extrapolated over-density for 
spherical collapse: it is 1.686 in the Einstein-de Sitter 
case, while it slightly d epends on redshift for m ore gen- 
eral cosmologies (e.g., iKitavama fc Sutd (|1996i )). The 
redshift Zf is the formation redshift of the halo, which 
for high mass halos is very close to the observed red- 
shift Zq. Hereafter we will thus make the approximation 
~ Zo = z. The second equality can be understood 
if we think of the density fluctuation as being time- 
independent while giving a redshift dependence to the 
collapse threshold Sc{zf) = Ac{zf)/ D{zf). Here D{z) 
denotes the general expression for the linear growth fac- 
tor, which depends on the background cosmology. In 
particular D{z) = (1 -f z)~^g(z)/g{0) where g(z) is the 
growth suppression factor for non Einstein-de Sitter Uni- 
verses. 

For Gaussian initial c ondit i ons, one obtains ( Kaiseil 
(fT98l : IPoirtzer fc Wisd (fl98l ): ljensen fc Szalavl (119861? ) 
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^h,M{r) = exp 
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-Ur) , (2) 



where an is the r.m.s. of the underlying dark matter 
fluctuation field smoothed on scale R. v ^ ^c/cfr and 



h,M{R) 



denotes the correlation function of the halos of 



mass M corresponding to radius R. Wr denotes the top- 
hat function of width R and the definition of <^ij(r) is 
J (Pr' £^{r')W^{\Y — r'l). In the second equality above we 
have expanded the exponential in series. The truncation 
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of the series holds for separations r ^ R. Thu s we obtain 
the well-known Kaiser's formula (jKaisedll 984( 1 of a scale- 
independent bias: 



(3) 



where 6z, = <^\- Compare d with the more refi ned re- 
lation s in e.g., iMo &: Whiti (|1996D and Catclan et all 
(jl998!) , an additive term l/(5c has been dropped as a 
consequence of the high-peak approximation in the first 
equality of Eq. 

Here the subscript L indicates that this should be con- 
sidered as a Lagrangian bias, because all correlations and 
peaks considered here are those of the initial density field 
(linearly extrapolated till the present time). Making the 
standard assumptions that halos move coherently with 
the underlying dark matter, one can obtain the final Eu- 
lerian b ias as = 1 + ^Li us i ng the techniques out - 
lined i n Efstathiou et aTl ri988fl. ICol e fc Kaiseil (|1989( ). 
IMo fc White (.19961 and lCatelan et al.. (,199^ 

The two-point correlation function of regions above 
a high threshold has been obtained, for t h e gen - 
eral non-Gaussian ca se, in jGrinstcin fc Wisd (fl98fih . 
iMatarrese et all (|1986D and lLucchin et all (119881 ): 
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j times (A^— j) times 



As before, for large separations we can expand 
the exponential to first order. T o leading order 
for non -Gaussianity of the tv pe (iSalopek fc Bondl 
1990t iGangu i et al.' [l99l : iVerde et all l200a 
Komatsu fc Spergel 200ll ) 



$ = + /nl - 



(5) 



where * denotes convolution, as in general /nl may be 
scale and configuration dependent, but for constant /nl 
it reduces to a simple multiplication. For simplicity, 
below we will carry out calculations assuming constant 
/nl and will generalize our results at the end. Here $ 
denotes Bardeen's gauge-invariant potential, which, on 
sub-Hubble scales reduces to the usual Newtonian pe- 
culiar gravitational potential, up to a minus sign. In 
the literature, there are two conventions for Eq. ([5]): the 
large-scale structure and the CMB one. Following the 
large-scale structure convention, here $ is linearly ex- 
trapolated at z = 0. In the CMB convention $ is in- 
stead primordial: thus /nl = g{z = oo)/ g{0)f§l^^ . In 
Eq. (O, (p denotes a Gaussian random field. For values 
of /nl consistent with observations, we can keep terms 
up to the three-point correlation function obtaining 
that the correction to the halo correlation function, A^^ 
due to a non-zero three-point function is given by: 



-44'^(xi,xi,X2) 



^j^^(xi,Xi,X2) 



(6) 



3. APPLICATION TO A LOCAL NON-GAUSSIAN 
MODEL 



2.0 
1.8 
1 .6 
1 .4 
1 .2 

1.0 

0.8 
0.0001 




0.0010 0.0100 
k 



0.1000 



1.0000 



Fig. 1. — The function Tuik) for three different masses: 1 X lO'^* 
Mq (sohd), 2 X 10^"' Mq (dotted), 1 X 10^^ Mq (dashed). 



We want to find an expression for the correlation func- 
tion of the late-time halos which form from the dark mat- 
ter over-density. 

In Fourier space, the present-time [z = 0) filtered lin- 
ear over-density 5r is related to $ by the Poisson equa- 
tion: 

Snik) = ^J^^WR(kMk) EE MnikMk) , (7) 



^ 'm,0 



where T{k) denotes the matter transfer function"' and 
VFfl(fc) is the Fourier transform of Wii{r). From Eq. ([5]) 
the definition of $ is 



m = 0(k) + /nl J ^0(k')0(k - k') , (8) 

whose bispectrum is 

B^ki, fc2, fca) = 2/nl [P,p{ki)P4k2) 

+ Pc^{k2)Pc^{k3) + P^{ki)P^{k3)] . (9) 

With these definitions the density bispectrum becomes 
Bs{ki,k2,k3) = MR{ki)MR{k2)MR{k3)B^{ki,k2,k3), 
where denotes the power-spectrum of the Gaussian 
field (f). The three-point function of Eq. ^ becomes 

e(3)(xi,Xi,X2)= (10) 
^/d3fcld3fc22/NLXfl(fcl)A^i?(fc2)A^fl(|ki+k2|)x 

[P4>{ki)P4>{k2) + 2P^(fci)P^(|ki + k2|)] e»(ki+k2)-(xi-^2) 
and Fourier transform of Eq. (fTU|) becomes 

^^Xfl(fc) j dkiklMR{ki)P^{ki) X 



d^lMR (^/^) [P^ (V^) + 2P^(fc)] 



where a = kl + P + 2kikfj.. 

^ The matter transfer function and the window functions cannot 
be neglected here: the initial conditions of Eq. JSJl are set out well 
before matter-radiation equality and the density field should be 
smoothed to define the halo mass. 
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Fig. 2. — The redshift dependence of Ahh/bh- 



4. RESULTS 



We can now write an expression for the non-Gaussian 
contribution to the halo power spectrum. Eq. ^ in 
Fourier space becomes: 

APh(fc) = hlj^^hi^5,P^s{k)TR{k) (11) 

where we have used 6o,l = ^clcf\, corresponding to the 
Lagrangian Unear bias that the halos would have in the 
Gaussian case, P^s{k) = A4R{k)Pcf,{k) and 



:FRik) = 



1 



dkiklMR{ki)P^{ki) x 



dfiM R ( Va) 



(12) 



The "form" factor J-R{k) is plotted as a function of k in 
Fig. [H for three different masses. 

The expression for the halo power-spectrum can be 
rewritten in a more convenient form where we can also 
make the redshift dependence explicit: 



Ph{k,z) - 



Sl{z)Pss{k,z) 



1 + 4/nl'5c(^) 



P^s{k)TR{k) 
Pssik) 



D\z)Ml{k)P^{k). 



where P55{k,z) = D^{z)Pss{k) — ±j v-iy^^i_RV 

We can now define the Lagrangian bias 6l of the halos 
from h\ — P\i{k,z)/ Pss{k,z) and use — 1 -I- 6^ to 
obtain the expression for the non-Gaussian halo bias 



Ac(^) 



1 + 2/nl 



A,(z) TR{k) 



D{z) MR{k) 



(13) 



Thus b{^^ = bh{l + Abh/bh) where bh denotes the 
halo bias for the Gaussian case. Abh/bh is 2/nl 
times a redshift-dependent factor Ac{z)/ D{z), plot- 
ted in Fig. times a k and mass dependent factor 
TR{k)/MR{k), shown in Fig. El 

5. DISCUSSION AND CONCLUSIONS 

We have obtained an analytic expression for the bias 
of dark matter halos for non-Gaussian initial conditions. 
The only approximations used in our approach are: i) 
high peaks, i.e. large values of v — 5c{z)/aR, (as in the 
original Kaiser's formula), which essentially amounts to 
a limitation on the mass range over which one can apply 



Fig. 3. — The scale dependence of Ab^/bh for three different 
masses: 1 X 10^"' Mq (sohd), 2 X 10^"' M© (dotted), 1 X 10^^ 
(dashed). 



this formula, and ii) large separation among the halos, 
which is the standard assumption allowing to use linear 
bias. While it is true that on large scales {k 0) the 
form factor, the transfer function and the window func- 
tion go to unity, on the scales of interest neglecting these 
terms may lead to errors on Abh and therefore on /nl of 
the order of 100%. Comparison of these analytical find- 
ings with simulations will be presented elsewhere (Grossi 
et al, in preparation). 

An advantage of our approach is that it can be eas- 
ily generalized to non-local and scale-dependent non- 
Gaussian models in which B^{ki, k2, fcs) is the dominant 
higher-order correlation and has a general form, obtain- 
ing 



Abh _ A,{z) 1 
bh D{z) STT^al 



dkik\MR(ki) X 
B^{ki,^/a, k) 



(14) 



Modeling the clustering of hot and cold CMB spots for 
non-Gaussian initial conditions, is a straightforward ex- 
tension of this calculation (Heavens et al., in prepara- 
tion). 

We envision that this calculation will be useful for con- 
straining non-Gaussianity from future surveys which will 
provide a large sample of galaxy clusters over a volume 
comparable to the horizon size (e.g., DES, PanSTARRS, 
PAU, LSST, DUNE, ADEPT, SPACE, DUO) or mass- 
selected (via the Sunyaev-Zel'dovich effect) large clusters 
samples (e.g., ACT, SPT). 
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